home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
Mac Magazin/MacEasy 19
/
Mac Magazin and MacEasy Magazine CD - Issue 19.iso
/
Wissenschaft & Technik
/
FFTs for RISC 1.1
/
ffttest.m
< prev
next >
Wrap
Text File
|
1996-01-30
|
342b
|
19 lines
N=64;
j=sqrt(-1);
a=ones(N,1);
row=1;
for i1=1:N;
a(i1) = sqrt(row+i1-1+.77777)+j*((row+i1-1+.22222)*(row+i1-1+.22222) / N - N/2);
end;
a(1) = N+3+j*imag(a(1));
a(2) = 1-N+j*imag(a(2));
A=fft(a);
a1=ifft(A);
if N<=32
A
a1
end
max(abs(A-atrans)./(max(abs(A))/N))
max(abs(a-aitrans)./(max(abs(a))/N))
%max(abs(a1-aitrans)./(max(abs(a1))/N))